Giant current-driven domain wall mobility in (Ga,Mn)As 
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We study theoretically hole current-driven domain wall dynamics in (Ga,Mn)As. We show that 
the spin-orbit coupling causes significant hole reflection at the domain wall, even in the adiabatic 
limit when the wall is much thicker than the Fermi wavelength, resulting in spin accumulation and 
mistracking between current-carrying spins and the domain wall magnetization. This increases the 
out-of-plane non-adiabatic spin transfer torque and consequently the current-driven domain wall 
mobility by three to four orders of magnitude. Trends and magnitude of the calculated domain wall 
current mobilities agree with experimental findings. 
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Spin-polarized currents can reverse the magnetization, 
excite spin-waves or move domain walls in ferromag- 
nets. These are intriguing phenomena which can be- 
come useful for magnetic memories. Ferromagnetic semi- 
conductors are especially interesting because the critical 
current density (j c ) required to move domain walls is 
two to three orders of magnitude smaller than in fer- 
romagnetic metals [2, 0, JjJfl- Much effort has been 

and theoreti- 
current-driven do- 
main wall dynamics [16j |. However, it is not understood 
why the critical current density for domain wall motion 
is so small in semiconductors compared to metals. 

The spin-current induced torque on domain walls can 
be written as a sum of an in-plane part and an out-of- 
plane part, T = T m + T ± . The in-plane (out-of-plane) 
torque is spanned by (perpendicular to) the gradient of 
the local magnetization. In ferromagnetic metals, the 
domain wall width (X w ) is large compared to the Fermi 
wavelength (Xf) and the spin-orbit coupling is weak. 
Here, the domain wall does not reflect electrons which 
adiabatically align their spin close to the local magneti- 
zation direction as they traverse the domain wall. The as- 
sociated angular momentum transfer induces an in-plane 

torque on the domain wall @, ©, m, m h eh, m 



-j s m x [m x (j ■ V)m]. 
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r m is also called the adiabatic spin torque. In Eq.([T]), j s 
is the spin current density. The unit vectors j and m are 
in the directions of the current and the local magnetiza- 
tion, respectively. At low current densities, t™ does not 
cause any steady-state motion of the domain wall IH, U | • 
However, even a small out-of-plane t orq ue t 1 - can induce 
a finite domain wall drift velocity |l2l . [T3L 14 1 . An of- 
ten used approximate form for the out-of-plane torque 
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-/3j s m x (j- V)m, 
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where (3 is a small dimensionless parameter. The domain 
wall drift velocity is controlled by the out-of-plane torque 
and is proportional to (3 [l2l.[l3j]. In literature, deviations 
from the adiabatic in-plane torque, e.g. t 1 - , are often 
called the non-adiabatic torque. We prefer to classify the 



torques as in-plane and out-of-plane, where the latter can 
induce a steady state domain wall motion. 

In this Letter, we show that the intrinsic spin-orbit 
coupling in the valence band of magnetic III-V semi- 
conductors strongly enhances the current-induced out- 
of-plane torque and thereby the domain wall velocity. 
t 1 - depends non-locally on the whole domain wall profile 
and cannot be described by the local expression given 
by Eq.((2]). We find that the steady-state domain wall 
velocity (v w ) is proportional to the current density (j). 
The relevant measure of current-driven domain wall mo- 
tion is the domain wall current-mobility /ij = v w / j. Us- 
ing realistic values for domain wall width and spin-orbit 
coupling, we find that /ij is enhanced by three to four 
orders of magnitude compared to a system with weak 
or vanishing spin-orbit coupling. This may explain an 
open question: why fij (or j c ) is much larger (smaller) in 
(Ga,Mn)As than in ferromagnetic metals [H, @, H, 0]- 

In conventional ferromagnetic metals only a very small 
fraction of the electrons contribute to the out-of-plane 
torque via spin-flip scattering at extrinsic impurities or 
via non-adiabatic corrections due to a finite domain wall 
width [TH, EH, 14 1 . Radically different from this, the 



strong intrinsic hole spin-orbit coupling in (Ga,Mn)As 
causes a finite domain wall resistance even in the adia- 
batic limit (X w 3> Xf), by preventing a large fraction of 
holes to adiabatically adjust their spins to the magneti- 
zation of the domain wall [3, E||. The intrinsic spin- 
orbit coupling in combination with the magnetization 
cause an anisotropic momentum distribution of propa- 
gating modes. Fig[T]a), b) and c) show distributions of 
transverse propagating modes in (Ga,Mn)As for differ- 
ent magnetization directions [3| ■ Only transverse modes 
that exist for all magnetization directions within the do- 
main wall, FigfT]d), will conduct. The remaining modes 
are totally reflected [3, Ei|. The fraction of reflected 
holes is large and comparable to the total number of con- 
ducting particles. This induces an enhanced domain wall 
resistance, as well as spin accumulation and mistracking 
between hole spins and the magnetization of the domain 
wall. We find that the total out-of-plane torque increases 
by three to four orders of magnitude compared to systems 
with vanishing spin-orbit coupling. The current-induced 
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mobility is subsequently enhanced approximately propor- 
tional to the total out-of-plane torque. 
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FIG. 1: a), b) and c) Distributions of transverse propagating 
modes in (Ga,Mn)As for different magnetization directions 
shown by the arrows. White (gray) color indicates that two 
(one) spin channels are open for transport, d) Distribution 
of transverse conducting channels through an adiabatic Bloch 
domain wall in (Ga,Mn)As. Arrows show the magnetization 
directions itinerant holes encounter while traversing a Bloch 
wall. For an adiabatic domain wall, transverse modes in a) 
that do not exist in d) are total reflected by the wall. Such 
reflections increase the out-of-plane torque on the domain wall 
which subsequently increases the domain wall drift velocity. 
The plots are created using Eq.© with 71=6.8, 72 = 2.7, 
h ex /e =l.5, e F /e =2.25. 

We shall use the simplest model that captures the es- 
sential physics of holes in (Ga,Mn)As, namely the 4-band 
Kohn-Luttinger Hamiltonian in the spherical approxima- 
tion 0, m ED, El 
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h(r) • S, (3) 



where, m and p are the bare electron mass and the hole 
momentum operator, respectively. S is a vector of 4 x 4 
dimensionless angular momentum operators for a 5=3/2 
spin, and 71 and 72 are Luttinger parameters. In Eq. © , 
h represents the exchange field from the localized mag- 
netic moments. The exchange field is assumed to have a 
constant modulus and to be homogeneous in the trans- 
verse directions, \h(y)\ = h ex . The parameter 72 deter- 
mines the strength of the effective spin-orbit coupling of 
the holes. For a given doping density, Eq.© with h=0 
describes hole doped GaAs and provides a Fermi energy 
e , a Fermi wavevector k —y/2m e a /h and a Fermi wave 
length A =27r/fc for heavy holes. Here, m =m/2.6 rep- 
resents the heavy hole mass in GaAs. 

Focusing on low temperatures and low bias transport 
properties, we consider the linear response at a Fermi 
energy e F . Our system is a discrete, rectangular conduc- 
tor of lengths L x , L y , L z and lattice constants a x ,a y ,a z 
sandwiched between two reservoirs. Born-von Karman 
boundary condi-tions are assumed in the transverse di- 
rections x and z. 

We model the dynamics of domain walls by the dimen- 
sionless Landau-Lifshitz-Gilbert (LLG) equation 



dm 



m x h 



dm 
~dT 



where m(y) should be understood as a unit macrospin 
for the transverse slice located at y, and a is the Gilbert 
damping constant. The dimensionless time is t = t/t 
where t =M s (\ /2ir) 2 /27A Here, M s , A and j=e/mc 
are the saturation magnetization, spin stiffness and the 
gyromagnetic ratio, respectively. The effective magnetic 
field reads 



h eff = V 2 m + h an ■ m - h pdXs ■ 



(5) 



where V=(A /27r)V. h an is the dimensionless anisotropy 
field, which is used to control the type and the width of 
the domain wall. Demagnetization fields for simple ge- 
ometries can be included in h an . The last term in Eq.© 
is the contribution to the spin-transfer torque from the 
non-equilibrium itinerant holes lH, 17 1. The dimension- 
less coupling constant h p d and the spin-density response 
function Xs are defined below. The LLG equation is nu- 
merically integrated using a fourth-order Runge-Kutta 
method [H. 

We calculate the current carryin g w ave function by 



(4) 



a stable transfer matrix formalism [18|, |23j. The non- 
equilibrium spin density (pS) ne is determined as a trace 
over all hole states between e F and e F + eV in the left 
reservoir. We define the linear response function for 
the spin density as % s = (e?ifc /jm )(pS)™ e , which gives 
h P 4=hh ex \ j /87rej4e , that is proportional to the current 
density j. Similarly, we also compute the non-equilibrium 
spin current density, j s , and the average spin per conduct- 
ing holes, (S(y))=}„e/jh. 

The parameters describing our conductor are a x /X = 
&z/X a = 0.2/27T, a y /X = 0.75/27T, L x = L z = 5la x and 
L y = 400a y . Furthermore, h ex /e = 1.5 and e F /e = 2.25. 
The ratio e F /h ex —1.5 is in the same order as experimental 
values [22j ]. The anisotropy/demagnetization field h x n =0, 
hy U =—l for simplicity and to ensure that the considered 
domain walls are Bloch walls, consistent with recent ex- 
periments 0, 0, We use h z n to control the domain 

wall width X w /X = ^ l/h% n /2ir. The Gilbert damping 
constant is assumed to be a=0.03 [2o| . 

The transfer of momentum and angular momentum 
from holes to the domain wall must be treated on an 
equal footing. We achieve this as follows: First, we find 
the equilibrium (j—0) domain wall configuration by in- 
tegrating Eq.© with h p d=0 to a sufficiently large time. 
The resulting equilibrium domain wall configuration is 
m x (y)=l/ cosh(y/X w ), m y (y)=0 and m z (y)=ta,nh(y/X w ), 
as expected from an analytical equilibrium solution of 
Eq.©. Next, a constant current is applied by choosing 
a finite h p d, and the following two steps are iterated: 1) 
compute Xs f° r a gi ven h=h ex m using Eq.©, and 2) 
integrate m = h/ h ex a time At using Eq.©. We use 
Ai=0.2, which is sufficient for convergence. 

Let us first consider a system without spin-orbit cou- 
pling, 72=0. We consider a X w /X =2 domain wall, close 
to the experimental values [l], HJ ■ Numerical calculation 



3 



(Sx) 


h z 






\T\ 

a) 

72 " - 






\ h z 








\ \T\ 

b) 

72=2 _ 


(St) J 


-15 -10 -5 


5 10 15 



0.01 



0.001 



1e-04 



1e-05 




2//Ao 



-20 -10 10 

2/Ao 



FIG. 2: Steady-state snapshots, a) and b): exchange field 
(h x /e and h z /e ), spin per conducting hole ({S x ) and (5 2 )) 
and absolute values of spin-torque per conducting hole \T\ as 
functions of position, for a) vanishing and b) finite spin-orbit 
coupling, c): absolute values of the out-of-plane torque \T y \ 
versus position for 72=0 (thick solid line), 72=0.5 (dashed 
line) and 7 = 2 (symbol +). For all plots, h p d= 0.001 and 
Au, / \q=2. 



of the conductance shows that only 0.1% of the incom- 
ing holes are reflected back by the domain wall. From 
FigI2k), we see that the hole spins follow the domain 
wall magnetization closely, and (S(y)) is virtually par- 
allel to h(y) throughout the system. Outside the wall, 
(S' 2 )«±1.25, since both the heavy (S z ~— 1.5) and light 
(S z ~— 0.5) holes participate in the transport. The abso- 
lute values of the dimensionless spin-torque per conduct- 
ing holes, T—(h/e ) x x s , is shown in Figl^k). It agrees 
well with the in-plane torque given by Eq.([T]) through 
the relation T=T m eA /2irj?i. The out-of-plane torque is 
very small, as expected for this case. 

As shown in FigfSb), turning on the spin-orbit coupling 
completely changes the physical picture. Numerical cal- 
culation of the conductance shows that the domain wall 
reflects 45% of the holes, causing spin accumulation and 
mistracking between carrier spins and the magnetization 
of the domain wall, particularly on the upstream side of 
the wall. Interference between incoming and reflected 
holes creates a spin-wave pattern in (S) and T, causing 
the reproducible "noise" in the figures. The shape of the 
domain wall is, however, virtually unchanged, due to the 
small current applied, h p< i — 0.001. From our previous 
study, we know that the majority of reflected holes con- 
sists of heavy holes [18(, in agreement with (S z ) « —3.0 
on the left side of the domain wall. 

As discussed in the introduction, the steady state do- 
main wall velocity is controlled by the out-of-plane torque 



12, l]| . Figf5t) shows the absolute value of the out-of- 
plane torque component, \T y \, for three different spin- 
orbit couplings. The spatial total of \T y \ denoted as 
(\T y \)— J dy\T y \ may be used as an estimate of the out- 
of-plane torque contributions to the domain wall drift 
velocity. Note that even when the spin-orbit coupling 
vanishes, the out-of-plane torque is finite since the sys- 
tem is not entirely adiabatic, A u ,/A =2 <C 00. Here, we 
find that \T y \ w 1% of the total torque |T|. When the 
spin-orbit coupling 72 increases from to 2.7, the spatial 
total out-of-plane torque (\T y \) increases by a factor of 
1000. This explains the dramatic increase in the domain 
wall drift velocity and thereby mobility. 
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FIG. 3: Domain wall displacement as a function of time for 
a) zero spin-orbit coupling for different domain wall widths, 
and b) fixed domain wall width for varying hole spin-orbit 
coupling strengths. For both plots, /i p( j=0.001. 

FigUM displays the domain wall displacement (x w ) as 
functions of time for varying domain wall widths for van- 
ishing spin-orbit coupling. x w completely saturates in 



the adiabatic limit, \ w / \ > 3, as expected 
Beyond the adiabatic limit, an increasing fraction of holes 
will be reflected, causing a out-of-plane torque on the do- 
main wall and a finite drift velocity, which, as can be seen 
in Fig[3£i), increases for decreasing domain wall width. 

Let us now consider a finite spin-orbit coupling. 
FigHJs) shows the domain wall displacement as a function 
of time for different spin-orbit couplings. We see that a 
larger spin-orbit coupling also increases the domain wall 
drift velocity due to the increased reflection of holes. The 
reflected holes increase the spin accumulation and mis- 
tracking, which, in turn, increases the out-of-plane spin- 
torque on the domain wall and thereby its drift velocity. 

FigHK) shows the domain wall drift velocity as a func- 
tion of h p d for varying spin-orbit coupling strength. We 
see that v w is proportional to h p d, i.e. the domain wall 
drift velocity is proportional to the current density, con- 
sistent with experimental observations [l[. In FigSJb) , we 
show the current-driven domain wall mobility as a func- 
tion of spin-orbit coupling. We see that the domain wall 
mobility increases four orders of magnitude when the hole 
spin-orbit coupling increases from 72=0 to a typical value 
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FIG. 4: a) Domain wall drift velocity as a function of h p d for 
varying spin-orbit couplings, 72. Note that the scale of the 
horizontal axis is different for each spin-orbit coupling, done 
to collect all graphs into a single plot, b) Current-driven 
domain wall mobility, /ij, and spatial total of the absolute 
values of the out-of-plane torque, (|7i/|), as functions of spin- 
orbit coupling. For both plots, \ w /\ =2. Lines are guides for 
the eye. 



of 72=2.7. In FigHb), we also show the spatial total out- 
of-plane torque as a function of the spin-orbit coupling. 
We see that (|^/(72)|)/(|^y(72=0)|) increases three orders 
of magnitude when 72 increases from to 2.7. Increasing 
spin-orbit coupling increases the out-of-plane torque and 
thereby the mobility. 

The effects of impurities [26| have been disregarded 
here. However, we believe that the present study still 
captures the essential physics. For example, we find that 
the intrinsic domain wall resistance persists in the dif- 
fusive transport regime 



27j. Therefore, we expect that 



the spin-orbit induced effects presented in this paper are 
important also for the current-driven domain wall dy- 
namics in the diffusive transport regime when the mean 
free path is not much smaller than X w . The expected im- 
purity induced reduction of the domain wall drift velocity 
is consistent with recent experimental findings [H, 0, El- 
Experiments find a two to three orders of magnitude en- 
hancement of the drift velocity, which is smaller than our 
computed three to four order of magnitude enhancement 
in the ballistic regime. 

In summary, the intrinsic spin-orbit coupling in 
(Ga,Mn)As increases the out-of-plane torque on the do- 
main wall and thereby the domain wall current-mobility 
by three to four orders of magnitude when the spin-orbit 
coupling increases from 72 = to an experimentally at- 
tainable value of 72=2.7. 
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